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We study the low-energy electronic properties of a junction made of two crossed metallic carbon 
nanotubes of general chiralities. We derive a tight binding tunneling matrix element that couples low- 
energy states on the two tubes, which allows us to calculate the contact conductance of the junction. 
We find that the intrinsic asymmetries of the junction cause the forward and backward hopping 
probabilities from one tube to another to be different. This defines a zero-field Hall conductance 
for the junction, which we find to scale inversely with the junction contact conductance. Through 
a systematic study of the dependence of the junction conductance on different junction parameters, 
we find that the crossing angle is the dominant factor which determines the magnitude of the 
conductance. 
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I. INTRODUCTION 

Since their discovery almost two decades agoi, carbon 
nanotubes (CNTs) have been the subject of intense ex- 
perimental and theoretical research. CNTs come in dif- 
ferent geometries: single wall tubes, multiwall tubes, and 
single wall tube bundles (ropes). Different CNTs are 
characterized by their chirality, which is a measure of 
how the tube lattice is oriented with respect to the tube 
axis. 

In order to understand the bulk behavior of CNT sys- 
tems, the electronic properties of different CNT geome- 
tries have been studied^, with some focus on the mul- 
tiwall tubes'^"^. These studies addressed the issue of 
intertube transport and its dependence on the chirali- 
ties of the tubes and the geometrical details. Most re- 
cently, films of CNTs have been used as a transparent 
conducting electrodeiS. Such films are made of networks 
of CNT ropes. In a rope, the constituent tubes have a 
distribution of chiralities. The resistance of these films 
is dominated by the tube-tube contacts(junctions). At 
such junctions, tubes with different chiralities intersect 
with different crossing angles. Therefore, it is desirable to 
study how these two factors (tube chirality and crossing 
angle) affect the transport properties of these junctions. 

Junctions of crossed CNTs have been studied 
experimentallyiiii^. In their work, Fuhrer et al^ii ex- 
plored the different electrical properties of different com- 
binations of crossed metallic and semiconducting tubes. 
Metal-metal junctions showed a contact conductance of 
the order of 0.02(4e^/ft,), despite the small junction area. 
This was attributed to the elastic deformation of the tube 
arising from the interaction with the substrate, result- 
ing in an intertube distance smaller than the expected 
3.4A. This increases the coupling at the contact re- 
gion, and provides a natural explanation of the relatively 



high measured conductance. Crossed Nanotube junc- 
tions have been theoretically studied for high symmetry 
cases onlyi^"— , where it was found that maximal contact 
conductance occurs when the tube lattices are commen- 
surate. Despite such effort, the general problem of two 
crossed CNTs has not yet been considered. Of interest is 
the effect of tube chiralities, the crossing angle as well as 
the relative orientation of the tube lattices on the trans- 
port properties of these junctions. 

In this paper we present a general study of that prob- 
lem in a tight binding framework. We derive a tunnel- 
ing matrix element that couples low energy propagating 
states on each tube. This matrix element depends on 
the chiral angles, the crossing angle, and the Fermi level 
of the junction. This allows us to systematically study 
the effects of the different junction parameters on the 
junction conductances. We find that due to the intrin- 
sic asymmetries of the junction, forward and backward 
tunneling between one tube and the other are generally 
unequal. Therefore, passing a current in one tube leads to 
the development of a non-zero voltage across the other 
one, thereby defining a zero-field Hall-like conductance 
for the junction. Furthermore, we find that this zero- 
field Hall conductance relates simply to the contact con- 
ductance of the junction. Our study also shows that the 
electronic properties of the junction sensitively-depend 
on the degree of matching between the tube lattices. For 
a given junction, this matching is controlled by the cross- 
ing angle, which results in an intertube conductance that 
varies by an order of magnitude for different angles. 

The paper is organized as follows. In section II we 
introduce our tunneling model and derive the tunneling 
matrix element coupling low energy tube states. In sec- 
tion HI we derive a formula for the junction conductance 
in terms of the microscopic conductances of the junction 
in a Landauer-Biittiker framework, and present some nu- 
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FIG. 1: Two crossed metallic tubes with different chiralities, 
and a crossing angle /3. 




FIG. 2: Details of the contact region. The dotted lines define 
the axes of the crossed tubes, with chiral angles 9i and 62- P 
is the crossing angle. Tubes face each other from the outside. 
The mismatch between the two lattices is parametrized by an 
angle (f) = 9i + 62 + P- The 2-D vector p defines the relative 
orientation of the origins of the two lattices(which are taken 
to be hexagon centers). 

merical results showing the dependence of the junction 
conductances on various junction parameters. The paper 
is concluded in section IV. 



II. TUNNELING MODEL 

In this section we describe the low-energy coupling be- 
tween two crossed tubes (Fig. [1} within a tight binding 
formalism. Nanotubes are assumed to be long and free 
of defects. The overlap region is shown in Fig. [2l where 
/? is the crossing angle, p is a vector defining the relative 
displacement of the origins(defined as hexagon centers) 
of the two lattices, and di and 62 are the chiral angles. 
Matching between the two lattices can by quantified by 
a registry angle 4> defined by 

= 0i + 02 + /3 (1) 

The uncoupled tubes are described by a nearest- 
neighbor tight-binding Hamiltonian 

^0 = - E E ^-''^a^CaJ (2) 
a=l,2 / . A 



where the index a labels the CNTs and (ij) is a sum 
over nearest-neighbor atoms on each CNT. The eigen- 
states of are plane waves localized on each CNT. The 
interaction Hamiltonian T-Lt is built in such a way that 
an electron can hop from any atom on one CNT to any 
atom on the other 

'Ht ^^Ujc\^C2j, (3) 

where tij depends on the positions and relative orienta- 
tion of the TT orbitals on atoms i and j, and varies ex- 
ponentially with the distance between the two hopping 
sites 

t,, -toe-''-/"", (4) 

where is a free parameter to be determined, and 
ao = 0.529A is the range of the tt orbitals. The inter- 
site distance dij is given by 

4 ^ {b + 2R- Rcos^ - Bcos^f 
^ R R 

+ (z2 sin /3 + i? sin i?sin — cos/3)^ 

^ R R ^' 

+ (z2 cos^ — 2i — i? sin ^ sin/3)^, (5) 

R 

where Za{ya) is the distance along the length (waist) of 
CNT a, 6 is the intertube separation at the point of clos- 
est contact. The CNTs are assumed to be of approxi- 
mately equal radii, which is denoted by R. We will also 
assume that hopping between the CNTs is dominated by 
Va, Za b, R. The total Hamiltonian of the system 

H = no + nT, (6) 

is then expressed in a plane wave basis by the transfor- 
mation 

^ai ~ ^fj^ ^ ! ^ ° '"C(j^(i)i(;„, (7) 

where Vai is the position vector of site i on tube a, rj spec- 
ifies the A 01 B sub-lattice, N is the number of graphene 
unit cells in the tube, and ka = {kay^ kax)- In this basis 

■Ho = -U ^ X!^k„c|,^k„CaBk„ +H.C., (8) 

a=l,2 ka 

where 

i=i 

and Aaj are the three nearest-neighbor vectors connecting 
the two sub-lattices of tube a. 

The Brillouin zones of the graphene sheets forming the 
two crossed tubes are shown in Fig. The wrapping of 
a graphene sheet corresponds to slicing the 2-D zone in a 
direction parallel to the tube axis. This yields a series of 
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FIG. 3: Brillouin zones of the lattices of the two crossed 
tubes. A. The zones are rotated with respect to each other 
due to the difference in chiral angles of the tubes, as well 
as the crossing angle /?. The axial momentum directions of 
the two tubes are shown. The rotation of one tube lattice 
with respect to the other is quantified by the angle <j). B The 
zones of the two lattices at one of the angles /3~ where the 
tunneling is dominated by pairs of Fermi points with zero axial 
momentum mismatch is zero and finite azimuthal momentum 
one. 



1-D bands, and in the case of a metallic tube, two bands 
cross at the Fermi energy. At low energy, one focuses 
on the neighborhood of the Fermi points (corners of the 
Brillouin zone), where k = dK.p + q, q = (0, g), a = ±1, 
andp = — 1, 0, 1 denote the three equivalent Fermi points. 
In this limit, the tube Hamiltonian is diagonalized by the 
transformation 



U ^ g-i(l/2)ae(T^g-i(7r/4)QO 



(10) 



where tr*''^ are the Pauli matrices, and we get the two 
eigenmodes ij^B^q (right movers, R) and i/'Lg (left movers, 
L). In this R/L basis 'Ho becomes 



EE' 

a=l,2 qa 



'"Pqai^lRg^tpaRq^ -^iLq^^aLqJ, (H) 



with vp — MTrd/2 and d is the nearest neighbor distance. 
The tunneling Hamiltonian "Ht can be expressed in the 
propagating states basis by following a similar procedure. 
Introducing the plane wave basis in Ht gives 



-Hi 



,iGi-(pi+»;iri)-G2-(p2+>)2T2) 



-EE- 

rjiVi GiG2kik2 
X iki+Gi.k2+G24iki'^';2k2 + ^-c- 



(12) 



where — a^Koi + q^, Pi are 2-D vectors defining the 
origins of each graphene lattice from the point of closest 
contact, r]i are ±1, are vectors joining the origin of 
each lattice to the nearest A sub-lattice site, and iki.k2 is 
determined from the 2-D Fourier transform of the spatial 
hopping element tij, and is given by 



*ki,k2 



t:, 



^{b+2R)fi^aoRf2 

L|sin/3| ' 



(13) 



where 



"■12 



^2z 



'2kizk2z cos/3 



(14) 



kly{k2z - kiz COS/S) + fc2y(fcl2 - k2z cos/3) 



+ 7;{{kl+kl)~{kl + kl)) 



sin /3 

ly + k2yJJ, 



tj 



27ra§6 



toe 



(15) 



(16) 



A is the area of a graphene unit cell, L is the length 
of the tubes, and tj is an overall magnitude of the k- 
space tunneling matrix element. This matrix element 
reproduces the result obtained for the case of parallel 
tubes(in the hmit (3 -> 0)^''. Taking the limits /3 ^ 
and i? — > oo gives us the hopping matrix element between 
two graphene sheets. This result can be used to calculate 
the low energy spectrum of two graphene sheets with 
Bernal stacking, as well as the transverse bandwidth of 
graphite^^. Fitting tj to the former, and for b = 3.4A, 
we get tj = 0.35 eV . 

We then express the tunneling Hamiltonian in the R/L 
moving basis. To do this, we notice from Eqs. (|13m5p 



that iki.k2 depends on the magnitudes of the momenta 
of the initial and final states, and decays rapidly as |k| 
increases. Therefore the sum over reciprocal lattice vec- 
tors in Eq. [12] is dominated by the first star. Ht thus 
becomes 

'Ht ^ ^ T{aiaiqi\a2(J2q2)lpla,aiqi'^2a2a2q2+H-C. 

(17) 



aia2cricr2Ql*?2 



where a = R,L, and T{aiaiqi\a2(72q2) is given by 



'■2£2 'P2 



T(aiaiqi|a2a2q2) = e'"^^"i -"^"^"^^^ 

ili2 = -l 

X ^QiKif +qi,Q2K2<>2 +q2-^'^0'iO-2 (18) 



where 



M 



fll ft2* 

J a\J — Q2 

ll fl2* 

a2 



fll fl.2* 

J a\J a2 
J ~a\J a.2 J—aiJ- 



ae 



and 







27r£ 
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(19) 



(20) 



(21) 



The form of tki,k2 (Eq- [T3lfT5|) shows that the cou- 
pling between propagating modes on the tubes depends 
on the geometry of the junction. In addition, the Fermi 
energy of the system has an effect on the coupling, as the 
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momentum difference between propagating states on tlie 
two tubes changes as the Fermi energy is changed. 

It should be noted here that our developed theory is a 
low energy one, and it applies where the linearization of 
the tube band structure is valid. An upper limit of the 
energy range (around the Dirac point) where our model 
is applicable is determined by the diameters of the tubes 
forming the junction. As an example, for tubes with 
diameters of ~ 1.4 nm, the higher sub-bands enter the 
picture at ~ 0.8 eV above and below the Dirac point^^. 
This energy scale decreases with increasing tube diam- 
eter (^ 0.3 eV for a diameter of ^ 4 nm). The low 
energy limit for the applicability of our model is imposed 
by the curvature gaps of small diameter tubes (^ 0.5 
nm), where the 7r-electron tight binding description fails 
to describe the band structure correctly close to the Dirac 



pointi 



.18.19 



These gaps are in the range of 50 meV for a 



diameter of about 0.7 nm^Si^. 

With Eqs. (|13ll2ip . we have a model that describes 
electronic coupling between two crossed metallic tubes. 
The virtue of our results is that they can be used to 
study the intertube conductances as a function of tube 
chiralities, crossing angle, Fermi energy, as well as the 
relative orientation of the tube lattices. Such a study is 
presented in the next section. 



III. JUNCTION CONDUCTANCES 

In this section, we derive a few formulas for the junc- 
tion conductances. A thorough treatment of the funda- 
mentals of mesoscopic transport is given by Dattap2^. 

One can abstractly view the two crossed tubes as a 
four-terminal device (FigH]). Coherent transport in such 
devices has been studied before^'^. In a four-terminal 
system with time reversal invariance, three different re- 
sistance measurements can be made. In Fig. [TJ one 
can imagine passing a current between terminals a and 
d, and measuring the voltage across terminals c and b, 
thus defining a resistance Rad.cb- Similarly, another two 
resistances, Racdb and Rab.dc, can be measured. The 
three resistances are subject to the simple constraint 
Rad,cb + Rac,db + Rab,dc = The forms of these re- 
sistances depend on the relations between different mi- 
croscopic conductances of the system under study. Ac- 
cording to the geometry of the system, one of these 
three resistances can be thought of as a zero-field Hall 
resistance(i?ab.(ic)- Four-terminal semiconductor hetero- 
structures have been previously studied2ii2^. The Hall 
resistance measured in such systems is due to an asymme- 
try induced by impurities. The CNT junction we study 
here is fundamentally different from those since the asym- 
metry causing the zero-field Hall resistance is an intrin- 
sic property of the junction. We will now derive ana- 
lytic results for the junction resistances in the Landauer- 
Biittiker formalism, and using our results from the last 
section, we will study the dependence of these resistances 
on different junction parameters. 



The transmission matrix between propagating states 
of different tubes can be calculated using Eqs. (|13ll2ip 
through Fermi's Golden rule. This allows us to calculate 
different intertube conductances. For linear response, 
these conductances are: 



Go 



(aa|HT|aV)|'Go, (22) 



where a, a' — R,L, pp is the density of states at the 
Fermi energy, and Go = 2e^ /h. Because of time re- 
versal symmetry, Gll = Grr and Glr = Grl, and 
therefore there are only two different intertube conduc- 
tances, denoted by Gf and Gb, respectively. In addition, 
in the absence of scattering, we have the intrinsic conduc- 
tances of the tubes, Gi = 2Go, i = 1,2. As we will see, if 
the distance between the two tubes is not very different 
from that in a bundle, coupling between the tubes will 
be weak, or equivalently, G/, Gt ^ Gi, G2. 

As pointed out, three conductance measurements can 
be done for the crossed tubes device. G^ = ^acbd 
first contact conductance of the junction, and is deter- 
mined by passing a current between terminals a and c 
and measuring the voltage across the h and c?(see Fig. 
[T]). G^.^ — R~^i,c ^^'^ other contact conductance of the 
junction, and is measured the same way as G^ but with 
c and d interchanged. The third conductance of the sys- 
tem, Gh = Rabcd^ ^ Hall-like conductance, where a 
current is passed in one tube and a voltage is measured 
across the other one. We use the Biittiker formula to ob- 
tain the following expressions for the three conductances: 



Gi'" 2{Gf + Gb) 2{2Go-Gbj) 



and 



{Gf - Gb) 



Gh 2(2Go — G/)(2Go — Gb) 



(23) 



(24) 



The three conductances satisfy a simple constraint^, 
Gl~^ - Gl^~^ = G^^ Since we expect G/,f, < Go, 
we can approximate the above formulas by 



1 



1 



1 



1 



G. 



Gi'" 2Gf + Gb 4Go 8G2 



and 



1 

Gh 



{Gf — Gb) 



8Gl 



(25) 



(26) 



where now the two contact conductances are equal to 
leading order in Gf/Go and Gb/Gg, and we will denote 
them by Gc- The product G^-Gh then becomes (to leading 
order in Gf_b/Go) 



GrGn — 



(27) 



where ?7= {G f—Gb) / {G f+Gb)- Since G/ and Gf, are gen- 
erally of the same order, 77 will be of order unity. Being 
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FIG. 4: Low energy band structure of two metallic tubes of 
different chiralities in the neighborhood of one set of K-points. 
The X-axis is along the first tube. The first tube bands (solid) 
and the second tube ones (dashed) are shown as a function of 
the axial momentum of the first tube, kix , so as to be able to 
compare the axial momentum mismatch of states at the same 
energy E. This results in the bands of the second tube being 
weighed by cos/3. At an energy Eo, the momentum mismatch 
for scattering between right moving states, AkiiR, is different 
from that between right and left moving ones, AkuL- 

a ratio between the difference and sum of the intertube 
conductances, the proportionality parameter 77 is inde- 
pendent of the tunneling strength io- Therefore, in such 
a system the Hall conductance is of order of the reciprocal 
of the contact conductance. This result is not restricted 
to our present system, but rather, it is a general result 
for any four terminal device which possesses time rever- 
sal symmetry, with two of its conductances much smaller 
than the other two. 

The zero-field Hall voltage that develops across a tube 
upon passing current in the second one is an intrinsic 
property of the tube junction. It is a manifestation of the 
difference between forward and backward transmission 
probabilities. Because of the Gaussian form of the matrix 
element (jl3p . it turns out that (for most chiralities and 
crossing angles) tunneling will be dominated by one set 
of Fermi points, say, Kio(of first tube) and K2o(of second 
tube). In that case, the tunneling matrix elements are: 

Ti+Rq,\+Rq2) - 2e''^'°-P'~''^-°-PH+K,^,+^,^+K,„+^, 

X COS — COS— , (28) 

and 

X COS — sm— . (29 

2 2 ^ ' 

The discrepancy between T{+Rqi\ + Rq2) and 
T{-\-Rqi \ -\-Lq2) arises from the momentum mismatch be- 
tween the initial and final tube states. To understand this 
more, one must have a closer look at the tube bands near 
the Fermi level. In Fig. |3]we show the lowest lying bands 



of a metallic tube(solid), and the same-energy states 
of another metallic tube of a different chirality (dotted). 
The horizontal axis is the momentum along the axis of 
the first tube, ki^. The second tube states(dashed) are 
plotted as a function of kiz as well. Tubes with different 
chiralities have different band-crossing points. The posi- 
tion of the Fermi level defines the axial momenta of the 
left and right moving states on each tube. The tunneling 
matrix element depends on terms of the form: 

Afc = — fc2z cos /3, (30) 

which is the difference between the axial momentum of 
one tube and the projection of the axial momentum of 
the second tube on the first one. At an energy Eq, the 
quantity Afc is given by: 

^kRL^RR. = — (l±cos/3) + (Kio)z-(K2o).cos/3 (31) 

and therefore iki ,k2 will be different for forward and back- 
ward scattering. 

It is worthwhile comparing the geometrical effects in 
the crossed tubes case to our previous work of two parallel 
tubesi^. For parallel tubes, the tunneling matrix element 
imposes axial-momentum conservation due to the effec- 
tive translational symmetry of the system, and approxi- 
mate momentum conservation in the azimuthal direction 
through its Gaussian dependence on the difference be- 
tween the final and initial states azimuthal components. 
For the case of crossed tubes, we have approximate 2- 
D momentum conservation; maximum coupling occurs 
when the 2-D momentum mismatch is minimum, which 
corresponds to highest matching between the tube lat- 
tices. We can summarize the effect of the crossing angle 
on the conductance in three cases. First, a general case 
where there is a finite 2-D momentum mismatch between 
the tube states(Fig. ^^). The greater this mismatch, 
the more suppressed is the inter-tube conductance. The 
second case is when the crossing angle is such that the 
Brillouin zones of the two lattices overlap, i.e. when 
(f) = 7i-|,n = 0,1,2, which will occur at crossing angles 
(defining j3 in the clockwise direction) 

/3+=n| + (0i+02). (32) 

At such angles, the 2-D momentum mismatch between 
the initial and final tube states is minimal. The third 
case is when the momentum mismatch is zero for one 
component only (Fig. [3j3). This occurs at crossing angles 
given by 

Pn^n^ + {9, - 92). (33) 

In general, conductance peaks occurring at these angles 
will be of lesser magnitude than those occurring at 
Therefore, we expect to see three peaks in the angular 
dependence of the intertube conductance corresponding 
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FIG. 5: Dependence of Gc on the crossing angle /3 for three 
junctions; I (17,2)-(10,10), II (17,2)-(15,6)(dashed), and III 
(13,7)-(13,7)(dotted), for Eo — 0. Different peaks in each plot 
mark the angles where there is high registry between the two 
tube lattices. The two arrows mark the angles /3~ and 13'^ for 
junction II. 



to the three angles /3+ where the tube lattice commen- 
suration is maximal, with some fine structure arising at 
the angles /3~ . 

Another effect of the crossing angle is that it controls 
the effective contact area between the two lattices, which 
is minimum at /3 = tt/2. Therefore, one expects conduc- 
tance peaks occurring close to a crossing angle of 7r/2 to 
be of a smaller magnitude than those occurring farther 
away from 7r/2. 

To illustrate the dependence of the junction conduc- 
tances on the different junction parameters, we study the 
roles of the chiral angles, crossing angle, relative orien- 
tation of the tube lattices, and the Fermi level of the 
junction. We perform this study for many junctions, of 
which we choose to show only three. The first (I) is a 
(17,2)-(10,10) junction, the second (II) is a (17,2)-(15,6), 
and the third (II) is a (13,7)-(13,7) one. The first two 
represent the general case of two tubes with different chi- 
rality, while junction III represents a junction which may 
be prepared experimentally by AFM manipulation^^. We 
take the intertube separation b to be similar to that of 
tubes in a bundle, b = 3.4A. The effect of the substrate 
would be to decrease the intertube separation. This can 
be easily incorporated into our model by using an effec- 
tive separation 5e// < 3.4A. For the relative orienta- 
tion of the CNT lattices, we study two limiting cases. 
The first is when the point of closest contact between 
the CNTs is two hexagon centers(HH-orientation), and 
the second is when that point is between an atom and a 
hexagon center (AH-orientation). 

We begin by discussing the dependence of the con- 
tact conductance, Gc, on different parameters. Figure [5] 
shows the angular dependence for Gc for the three junc- 
tions. The Fermi energy, £^0 = 0, is at the K point, and 
the CNT lattices are oriented with an HH orientation. 
For junction I (solid plot), we see that there are three 
main peaks in the conductance, which can be related 
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FIG. 6: Dependence of Gc on the crossing angle /3 for 3 junc- 
tions. The radii of the tubes in the 3 junctions are R(solid), 
2R(dashed), and 3R(dotted), and R ~ 1.5 nm. The CNTs are 
assumed to have an HH-orientation contact, and the Fermi 
energy Eq = 0. 



to maxima in the transmission probabilities between the 
two tubes at these crossing angles. As argued before, 
and since the tubes are of different chirality, there exist 
certain crossing angles at which lattice commensuration 
occurs, and in such cases the transmission probabilities 
will be largest, leading to a conductance maximum. Ac- 
cording to Eq. [321 the crossing angles where we should 
see such peaks are /3~^ = 36°, 96°, and 156° for junction 
I, which is what we see in Fig. [Sj On the other hand, 
crossing angles which maximize the lattice mismatch will 
result in a contact conductance which is roughly an order 
of magnitude lower. Junction II shows a qualitatively- 
similar behavior(Fig [S] dashed plot). For this junction, 
lattice matching occurs at l3+ = 22°, 82°, and 142°. 
For junction III, /3+ = 40°, 100°, and 160°. We also 
see that for this junction the conductance increases as 
/3 — )■ 0°,180° as the system becomes periodic and the 
tunneling matrix element takes a (5-function formic. 

The discrepancy in the magnitudes of the three con- 
ductance peaks (for each junction) can be understood in 
terms of the contact area of the junction. The transmis- 
sion probability, and hence the conductance, increases 
with the area of overlap between the two lattices. This 
area is minimum at perpendicular crossing, which makes 
the conductance peaks closer to /3 = 0°,180° relatively 
larger than those near /3 = 90°. 

The slight deviation of the conductance peaks from 
the lattice-commensuration angles that we see in Fig. [5] 
is due to the finiteness of the tube radii. For a radius 
of ~ 1.5 nm, the curvature causes a deviation of a cou- 
ple of degrees. This deviation decreases with increasing 
tube radii. In addition, this deviation is smallest at per- 
pendicular crossing, as the curvature effects are minimal. 
We show that in Fig. [51 for junction II(dotted plot) and 
two other junctions of the same chirality but with bigger 
radii; a (34,4)- (30, 12) junction (radii ~ 3 nm, dashed), 
and a (51,6)-(45,18) junction (radii ~ 4.5 nm, solid). As 
the radii increase, deviation from the angles Z?"*" becomes 
negligibly small. We also see that the peaks occurring at 
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FIG. 7: Dependence of Gc on the crossing angle /3 for junction 
I(17,2)-(10,10) with HH-(solid) and AH-(dashed) orientations 
at Eo = OeV. 
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FIG. 9: Dependence of rj on the Fermi energy for junction 
I(17,2)-(10,10)(sohd) and junction III(13,7)-(13,7) (dashed) 
at /3 = 40° with HH orientation. 
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FIG. 8: Dependence of Gc of junction I(17,2)-(10,10) on the 
crossing angle /3 for different energies: Eq = OeV(solid), and 
Eo = -0.2e'l/(dashed). 



angles [3~ get increasingly small as the radii of the tubes 
increase. These are the angles at which there is a finite 
mismatch in only one of the 2D-momentum components, 
which causes these peaks to be suppressed as the tube 
radii are increased. 

The dependence of the conductance on the relative ori- 
entation of the tube lattices is investigated in Fig. [71 
where we show the angular dependence of the conduc- 
tance of junction I for two orientations; HH, and AH. A 
first observation is that the relative orientation has a neg- 
ligible effect on the positions of the conductance peaks(as 
should be the case since these positions are a function of 
the Fermi-points mismatch). In addition, it is clear that 
the HH orientation gives a higher conductance, specially 
at the peaks where the lattice commensuration occurs. 
This is expected as in the HH orientation more lattice 
sites are matched, compared to the AH one. 

Whereas varying the crossing angle of a certain junc- 
tion can be a challenging experimental task, changing 
the Fermi level should be easily achievable through elec- 
trostatic doping. Therefore, the Fermi level dependence 
of the contact conductance is of interest. In Fig. |51 we 
show the angular dependence of the contact conductance 
for junction I for different energies. We find that whereas 



changing the Fermi level has a negligible effect on the po- 
sition of the conductance peaks, it does affect the mag- 
nitude of the conductance. This can be understood in 
band structure terms. Tuning the Fermi level away from 
the K points changes the momentum mismatch between 
the initial and final states involved in the tunneling(Eq. 
I3ip . thereby changing the transmission probability and 
the conductance. Another quantity which can be probed 
experimentally is rj. The virtue of such measurement is 
that, being readily predictable from our theory, it pro- 
vides a direct test for it. As mentioned before, rj is inde- 
pendent of the intertube separation. Our model predicts 
r] to be of order unity. The variation of 77 with the Fermi 
level is also achievable experimentally. In Fig. |9l we 
show the variation of 77 with the Fermi level for junctions 
I(solid) and Ill(dashed), for /3 = 40°. As we see, rj is of 
order unity. We find this to be the case for most junc- 
tions. A sign change of rj reflects the fact that the Fermi 
level affects the forward and backward conductances dif- 
ferently. The unique point in parameter space that has 
77 = is a point of high symmetry in the sense that the 
many asymmetries of the junction counter-act to give 
equal forward and backward transmission probabilities. 

Although electron-electron interactions are not in- 
cluded in our treatment, Eq. 1^7^ still holds. Such 
interactions tend to renormalize the tunneling density 
of states^''; p^'* = Pf(^)", and therefore 77 would 
not change if such interactions were included. The es- 
timates a ~ 0.6 and Tp ~ leV give p™* ~ O.lpp 
at room temperature. Another effect that we have ig- 
nored is the small deformation occurring in the tubes at 
the contact region, as predicted by molecular dynamics 
simulations^-'-^. Such deformations slightly increase the 
contact area between the two tubes, and may also cause 
the on-tube conductance to be slightly lower than the as- 
sumed 2Go value due to possible back-scattering. These 
two effects are expected to be of more relevance for large 
diameter tubes, where faceting at the contact region is 
more pronounced. Detailed treatment of the deformation 
effects requires a microscopic approach using a real space 
Hamiltonian^. 
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IV. CONCLUSION 

In this paper, we study the low-energy electronic prop- 
erties of a junction formed by two crossed metallic CNTs 
in a tight-binding framework. We derive a tunneling ma- 
trix element that couples the low-energy electronic states 
on the two tubes. Tunneling is found to be approximately 
momentum conserving in the sense that it has a Gaus- 
sian dependence on the momenta of the initial and final 
states. The magnitude of the coupling is determined by 
the intertube separation, and the crossing angle (which 
determines the contact area of the jmiction). 

The developed model allows for a clear understanding 
of the effects of the different junction parameters on its 
low energy electronic properties. We find that the in- 
trinsic asymmetries of the junction create a discrepancy 
between the forward and backward hopping between the 
tubes. This defines a zero-field Hall-like conductance for 
this four-terminal device. Using a Landaucr-Biittiker for- 
malism, we calculate the different conductances of the 
four-terminal junction. We find that the contact con- 
ductance scales inversely with the zero-field Hall conduc- 
tance of the junction. 

We also find that, in general, the crossing-angle depen- 
dence of the junction contact conductance Gc exhibits 
three peaks over the angular range. These peaks corre- 
spond to angles where matching between the tube lat- 
tices is highest, thereby maximizing the coupling. The 
relative magnitude of these peaks is understood in terms 



of the effective contact area between the tubes. Rela- 
tive orientation of the tube lattices is found to affect the 
strength of the coupling between the tubes, though its 
effect on the contact conductance is small compared to 
that of the crossing angle. We also find that the contact 
conductance varies with the Fermi level of the junction, 
but also in a way which is less dramatic than that of the 
crossing angle. We therefore conclude that the two cru- 
cial parameters in determining the contact conductance 
are the tube chiralities and the crossing angle. 

We believe that this work has important implications 
on studies involving CNT networks, where it is custom- 
ary to assume a fixed contact conductance between var- 
ious tubes in a network. With our model, and given a 
certain distribution of chiralities in the network, one can 
study the effect of the chirality distribution on the CNT 
network resistance. This may prove that some chirality 
distributions are favored over others. 
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